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Modules and packages for the one-loop calculations at partonic level represent the first level of SANC 
output computer product. The next level represents Monte Carlo integrator mcsanc, realizing fully differential 
hadron level calculations (convolution with PDF) for the HEP processes at LHC. In this paper we describe 
the implementation into the framework mcsanc first set of processes: DY NC, DY CC, fifi -> HW^{Z) and 
single top production. Both EW and QCD NLO corrections are taken into account. A comparison of SANC 
results with those existing in the world literature is given. 



1. INTRODUCTION 

Recent reviews of theoretical predictions and their 
uncertainties for basic LHC processes in the Standard 
Model (SM) can be found in Reports [T] and [3]. 

The interpretation of high-quality data of the LHC 
demands an equally high precision in the theoretical pre- 
dictions at the level of quantum corrections. Apart from 
a detailed knowledge of higher-order EW and QCD cor- 
rections, the combination of their effects must be inves- 
tigated. Advanced computational tools were developed 
to control the interplay of EW and QCD corrections: 

la, a, 0, m and h. 

In this paper new results of the computer system 
SANC (Support of Analytic and Numerical Calculations 
for experiments at Colliders) [8] are presented. In this 
system it is possible to achieve the one loop level predic- 
tions in the EW and QCD sectors on the same platform 
of the analytic procedures. 

The first level of the computer products SANC 
are: analytical modules for scalar Form Factors (FF) 
and Helicity Amplitudes (HA) and accompanying 
bremsstrahlung contributions (BR or MC) and the 
s2n.f package producing the FORTRAN codes [5]. 

In this paper we discuss in some detail the results 
of the implementation at hadronic level in the newly 
developed mcsanc-1.0 integrator, based on the above 
mentioned modules. The processes are marked by pro- 
cess identifiers: pid=cnn, c=charge: 0-NC, ±-CC, and: 
nn=01(e), 02(^), 03(t) etc, see below. 

• Drell-Yan-like single W production: pid = ±102. 



(1) 



• Drell-Yan-like single Z production: pid = 002. 

q + q^l+ + r (2) 

• HW^{Z) production: pid = ±104 (004). 
At the parton level we consider 



fJ[HW^{Z) ^ 



(3) 



(where /i stands for a massless fermion of the SM, while 
specifically for bosons we use Z, W^, H). It should be 
emphasized also that the notation ffHW — means 
that all external 4-momenta flow inwards; this is the 
standard SANC convention which allows to compute 
one-loop covariant amplitude (CA) and form factors 
(FF) only once and obtain CA for a specific channel 
by means of a crossing transformation. 

• the s and t channels of single top quark production: 
pid = ±105(s),±106(i). 



tbud — ^ and tbud — >■ 0. 



(4) 
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Previous studies of these processes by SANC system, 
i.e. creation of the analytic platform and modules at the 
parton level were presented in [8l ITO l [TT l [T2 l [T3 l [M l [T5 ] . 

The paper is organized as follows. First, an overview 
of the mcsanc-1.0 integrator is given in section 3. In 
subsection [21] we describe a list of contributions to hard 
sub-processes and introduce their enumeration. Then 
we describe the parallel calculations issue. Section 3 
contains numerical results for the processes in QCD and 
EW sectors. Input parameters, kinematical cuts, and 
the used PDFs can be found in subsection [3T] Fur- 
ther, we present the complete predictions for inclusive 
cross sections at LO and NLO levels in the EW and 
QCD sectors for processes (1I])-(|11). We systematically 
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compare our results for NLO QCD corrections with the 
program MCFM [16] , [17] and, whenever possible, with 
other codes existing in the literature: [TB], [T^ (EW) 
and [10] (QCD). 

In section 4 we summarize our results. 



2. SANC INTEGRATOR 

21 Description of id's for hard sub-processes 

At NLO level several hard sub-processes contribute 
to a given process. In general, it consists of sev- 
eral parts: LO-lowest order, Virt- virtual. Real-real 
brems(glue)strahlung and Subt-subtraction; Real, in 
turn, is subdivided into Soft and Hard contributions. 
We enumerate them through id=0-6: 



idO: LO, 2 
idl 



2, tree-level, qq' NC or CC sub-processes. 



Subt term, responsible for the subtraction of the 
initial quark mass (m^) singularities for qq' sub- 
processes, computed in a given subtraction scheme 
(MS or DIS). It depends on In(mq). 

id2: Virt represents only the NLO EW parts, stands 
for pure EW one-loop virtual contributions. It de- 
pends on niq and may depend on an infrared regu- 
lator (e.g. on an infinitesimal photon mass). It is 
not present for QCD NLO contributions, where it 
is added to the soft contribution (see next item). 
For DY NC NLO EW process it contains all vir- 
tual contributions, both EW and QED. 

id3: For all processes, except DY NC, this stands for 
the sum of Virt and Real Soft (QED/QCD) con- 
tributions, therefore it does not depend on the in- 
frared regulator but depends on ruq and on the 
soft-hard separator w. For DY NC NLO EW pro- 
cesses it is just the Real Soft QED contribution 
that depends on the infrared regulator, on and 
on the soft-hard separator lj. 

id4: For all processes this is just the Real Hard 
(QED/QCD) contribution that depends on rrig 
and on ui. 

id5: Subt term is responsible for the subtraction of 
the initial quark mass singularities for gq{gq') sub- 
processes (also computed in MS or DIS schemes). 
It contains logarithmic singularities in m^. 

id6: The gluon-induced sub-process-an analog of id4 
for gq{gq') sub-processes. They also contain loga- 
rithmic mass singularities which cancel those from 
id5. 



The quark mass is used to regularize the collinear 
divergences, the soft-hard separator is a remainder 
of infrared divergences. The sum of contributions 
with id3 and id4 is independent of a). The sums 
idl+id2+id3+id4 and id5+id6 are separately indepen- 
dent of mq. Therefore, the entire NLO sub-process is 
independent of both unphysical parameters w and mq. 



22 Parallel calculations 

The mcsanc program takes advantage of paralleliza- 
tion in the Cuba library [21] , [22] , used as a Monte Carlo 
integrating tool. However, the parallelization efhciency 
is reduced by the overhead of inter-process communica- 
tions. 

Figure [T] shows time required to complete the NLO 
EW cross section calculation depending on the num- 
ber of active CPU cores. The test was run on a dual- 
processor Intel ® Xeon ® machine with 12 real (24 vir- 
tual) cores with Linux operating system. The upper 
plot summarizes multicore CPU productivity: "total" 




12 16 20 24 

number of cores 

Figure 1. CPU usage and load efficiency for the mcsanc 
program depending on the number of processor cores 



is the wall clock time passed during the run; "user" is 
the CPU time consumed by the program (roughly equals 
to wall clock time multiplied by the number of cores in 
case of 100% efhciency); "system" is the time spent by 
the operating system on multiprocessing service. 

One can see that the parallelization is efficient with 
number of cores up to 8, after which the total run time 
does not significantly decrease and the overhead CPU 
time ("user") grows. It is also apparent from the lower 
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plot that the average CPU load efficiency drops below 
50% with more than 8 cores active. 



3. NUMERICAL RESULTS 

In this section the results, obtained by the mcsanc 
integrator, realizing fully differential hadron level calcu- 
lations for the processes ©-(HI) are presented. 

We produce comparison with numerical results for 
the NLO QCD corrections for all our processes between 
SANC and (TB]. For the NLO electroweak corrections 
for Drell-Yan NC and CC processes (pid = 002, ±102) 
this was done early within the workshop [23]. For WH 
production, i.e. pid = ±104, in EW sector between 
SANC and flU and in QCD sector between SANC and 
[20] for pid = ±105, ±106. 



31 Setup: PDF, cuts, input parameters 

For the numerical results in this section we have used 
the following setup. 

• PDF set, scales, as- We use CT10(f [scale] ) PDF 
from the LHAPDF library and compute Qs via a call 
to alphasPDF(r [scale] ) . Usually we set factorization 
scale (f scale) equal to renormalization scale (rscale) 
and different for the processes under consideration: My 
for DY-like single V production; My+H for the pro- 
cesses Eq.(j31); mt for the processes Eq.(|3]). 

• Phase-space cuts. We use loose cuts: for the fi- 
nal state particle transverse momenta pt > 0.1 GeV, 
no cuts for their rapidities and for the neutral current 
DY, in addition, AIi+i- > 20 GeV. We demonstrate nu- 
merical results for Drell-Yan only for muon case and 
we are not dealing with effects of recombination. We 
choose CO = 10"'' and the cms energy ^/sq = 14 TeV if 
not stated otherwise. 

• Set of EW scheme and input parameters. We 
choose the EW scheme, and input parameters are 
taken from PDG-2011 (on 16/05/2012): 

Coupling constants: a = 1/137.035999679, Gp = 
1.16637 X 10"^ Boson masses: Mw = 80.399 GeV, 
Mz = 91.1876 GeV, Mh = 120 GeV. Boson widths: 
Tz = 2.4952 GeV, = 2.085 GeV. CKM matrix: 

Vud = 0.9738, Vus = 0.2272, Krf = 0.2271, V^s = 0.9730. 
Lepton masses: rrie — 0.510998910 MeV, = 
0.105658367GeV, = 1.77682 GeV. Heavy quark 
masses: mi, — 4.67 GeV, mt = 172.9 GeV. Masses of the 
four light quarks are taken from [24]: = 0.066 GeV, 
m„ = 0.066 GeV, = 0.150 GeV, m^ = 1.2 GeV. 



32 Example of M^+^- distributions, DY NC 

The standard ATLAS Monte Carlo (MC) genera- 
tion uses the PYTHIA-PHOTOS chain of programs. 
PYTHIA [33], [H] has the leading order (LO) matrix 
element for a given process and takes into account Par- 
ton Showers (PS); PHOTOS [37], [33] and [H] describes 
multiphoton emission from the Final State (FSR) di- 
lepton system. This procedure does not include certain 
next-to-leading order (NLO) EW corrections, like Pure 
Weak (PW) contributions, Initial-Final QED interfer- 
ence (IFI) and what remains from Initial State Radi- 
ation (ISR) after subtraction of collinear divergences. 
In SANC one can evaluate the entire effect of these 
corrections as a difference of complete NLO EW cor- 
rections and the QED FSR corrections. See, for ex- 
ample, the distribution of the complete NLO EW cor- 
rections over Z boson invariant mass, J(M^+^-), for 
Drell-Yan-like single Z production around Z resonance, 
{S = {da^^° /dM^+^-)/{da^° /dM^+f,-) - 1 and with 
taking into account only QED FSR corrections. Fig- 
ure 2. 



PP ^ ^V(Y) (Vs^ = 7 TeV) 
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Figure 2. 5 in % with complete NLO EW (solid his- 
togram) and FSR (dashed histogram) corrections 
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Figure 3. Difference correction 5^^'^ in % for tire distri- 



bution over M^+^~ 
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NLO and FSR distributions in Figure 2 
arc barely distinguishable. The difference 

jDiff ^ jNLO-FSR(-^^^^_) jg gjjQ^jj Figure 3. 

As is seen, for the M^+^- interval around the Z res- 
onance, (5^'*^ varies from +5% at the lower edge to —1% 
at the upper edge and therefore cannot be neglected, if 
the precision tag is equal to 1%, say. 

33 Numerical results and comparison for LO, 
NLO EW, NLO QCD RC 

• In Tables 1-3 we present LO and NLO inclusive 
cross sections for processes (1-2), (3) and (4), respec- 
tively. 



pid 


002 


102 


-102 


LO 


3338(1) 


10696(1) 


7981(1) 


LO MCFM 


3338(1) 


10696(1) 


7981(1) 


NLO QCD 


3388(2) 


12263(4) 


9045(4) 


NLO MCFM 


3382(1) 


12260(1) 


9041(5) 


^QCD 


1.49(3) 


14.66(1) 


13.35(3) 


NLO EW 


3345(1) 


10564(1) 


7861(1) 


Sew 


0.22(1) 


-1.23(1) 


-1.49(1) 



Table 1: NC and CC DY processes, i.e. for pid 
= 002, ±102. LO, NLO EW, NLO QCD cross sections 
are given in picobarns and compared with correspond- 
ing values obtained with the aid of the program MCFM. 
Also correction factors are shown in %. The numbers il- 
lustrate good agreement within statistical errors of MC 
integration. 



pid 


004 


104 


-104 


LO 


0.8291(1) 


0.9277(1) 


0.5883(1) 


LO MCFM 


0.8292(1) 


0.9280(2) 


0.5885(1) 


NLO QCD 


0.9685(3) 


1.0897(3) 


0.6866(3) 


NLO MCFM 


0.9686(1) 


1.0901(2) 


0.6870(1) 


Sqcd 


16.81(3) 


17.47(3) 


16.72(5) 


NLO EW 


0.7877(1) 


0.8672(2) 


0.5508(1) 


Sew 


-5.00(2) 


-6.52(2) 


-6.38(3) 



Table 2: The same is in Table 1 but for processes of 
HZiW"^) production, i.e. pid= 004, ±104. 



• In Table 4 we show QCD and EW cross sec- 
tion contributions to the processes of iJW^ production, 
pid= ±104, detailed over id's of the mcsanc integrator. 
As is seen for the chosen setup (Subsection 3.1) there is 



pid 


±\JO 


-1 (Ti 

±\JO 


LO 


5.134(1) 


3 205CLI 


LO MCFM 


5.133(1) 


3.203(1) 


NLO QCD 


6.921(2) 


4.313(2) 


NLO MCFM 


6.923(2) 


4.309(1) 


6qcd 


34.79(5) 


34.56(8) 


NLO EW 


5.022(1) 


3.140(1) 


6ew 


-2.18(1) 


-2.02(2) 


pid 




1 rifi 

-lUU 








LO MCFM 


158.69(7) 


95.27(4) 


NLO QCD 


152.13(9) 


90.44(7) 


NLO MCFM 


152.07(14) 


90.50(8) 


Sqcd 


-4.17(6) 


-4.08(8) 


NLO EW 


164.44(5) 


98.65(4) 


Sew 


.^.59(3) 


3.()G(r)) 



Table 3: The same is in Table 1 but for single top, s 
and t channels, i.e. for pid= ±105, ±106. 

strong cancellation of the Soft and Hard contributions 
id's=3,4 and gluon induced contributions id's=5,6, 
the sum of contributions id's=5+6 being negative. We 
remind that the sum of all contributions is independent 
of the unphysical parameters u> and m,. 



QCD 


pid 


104 


-104 




idO 


0.9277(1) 


0.5883(1) 




idl 


0.6916(1) 


0.4860(1) 




id3 


-10.9233(1) 


-6.9139(1) 




id4 


10.4547(1) 


6.5737(1) 




id5 


-0.9717(1) 


-0.6733(1) 




id6 


0.9107(1) 


0.6258(1) 




NLO QCD 


1.0897(3) 


0.6866(3) 


EW 


pid 


104 


-104 




idO 


0.9277(1) 


0.5884(1) 




idl 


0.0100(1) 


0.0070(0) 




id2 


-0.0560(0) 


-0.0349(0) 




id3 


-0.1592(1) 


-0.1003(0) 




id4 


0.1448(1) 


0.0907(1) 




NLO EW 


0.8672(2) 


0.5508(1) 



Table 4: EW and QCD radiative corrections in pi- 
cobarns detailed over id's for pid= ±104, parameter 
w = 10-^ 
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• The other comparisons. 

1) EW corrections: a comparison of EW and QCD 
NLO corrections between mcsanc and papers |18) and 
[T5] was done using corresponding setup. We received 
good agreement within statistical errors with the results 
presented in the references. 

2) QCD corrections: a comparison between inclusive 
LO and NLO cross sections from mcsanc and Table 1 
of paper [^D] was carried out using tuned setup. Agree- 
ment within MC errors was found for LO, while for NLO 
only a qualitative agreement was reached, since we did 
not manage to reproduce the corresponding value for 
as{r) 

More comparisons, including differential distribu- 
tions (as in papers [4], [30]) will be presented elsewhere. 

4. CONCLUSIONS 

To match the experimental accuracy at the LHC, 
we direct our effort to developing a programming en- 
vironment for the calculation of processes at one loop 
level and to creating the mcsanc integrator with EW 
and QCD branches at hadron level. 

In this paper we have presented results for EW and 
QCD corrections to the following processes: Z and W 
production, HZ and HW^ production, and single top 
production processes. Our investigation confirms that 
NLO precision level and combination of EW and QCD 
corrections are mandatory for the precision tag of the 
LHC experiments. 

We thank A. Arbuzov, G. Nanava and A. Ochirov 
for essential contributions to the mcsanc project, and 
V. Kolesnikov and R. Sadykov for useful discussions. 

This work is partly supported by Russian Founda- 
tion for Basic Research grant N° 10-02-01030-a. A. 
Sapronov is cordially indebted to "Dinastiya" founda- 
tion 2011-2012. 
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